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Starting with two weakly-coupled anti-ferromagnetic spinor condensates, we show that by chang- 
ing the sign of the coefficient of the spin interaction, (72, via an optically-induced Feshbach resonance 
one can create an entangled state consisting of two anti-correlated ferromagnetic condensates. This 
state is maximally entangled and a generalization of the Bell state from two anti-correlated spin— 1/2 
particles to two anti-correlated spin— N/2 atomic samples, where N is the total number of atoms. 



Quantum entanglement is the single most fundamen- 
tal difference between quantum mechanics and classical 
mechanics [l]. It is also the key ingredient that enables in- 
formation transmission and processing beyond that pos- 
sible classically Q- Quantum degenerate gases 0,0, due 
to their purity, weak environmental coupling and high 
experimental control, are a natural candidate to achieve 
such entangled states. In particular, with the experi- 
mental success of all optical trapping of a Bose-Einstein 
condensate , it is now possible to explore atomic gases 
with spin degrees of freedom. There have already been 
a number of proposals to create multi-particle entangle- 
ment between atoms with an internal degree of freedom 
El H> IH • These proposals have concentrated on creat- 
ing spin-squeezed states 10] or pairwise entanglement 11], 
in contrast, in this letter we wish to propose a method to 
produce a maximally entangled state of the spin degree 
of freedom of two spatially separated spinor condensates 
via an anti-ferromagnetic to ferromagnetic transition. 

The original maximally entangled Bell state consists of 
a pair of anti-correlated spin— 1/2 particles 
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We can also straightforwardly generalize the above state 
to an arbitrary spin—/ particle. This generalized Bell 
state is also maximally entangled and has the form 
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such that each component of the spin, mi = —I, • ■ ■ I, 
of particle 1 is anti-correlated with that of particle 2. 
The "particles" we consider in this paper are in fact 
two ferromagnetic spin-1 atomic samples localized at 
the minimums of a double-well optical potential, where 
the spin variable I — N/2 is equal to the number of 
atoms in each well and can have the mesoscopic values 
~ 10 3 — 10 7 4 : ]. Measurement of the entangled state is 
performed in a similar way to the Stern-Gerlach experi- 
ment, where a magnetic field gradient applied to one of 
the atomic samples will spread the angular momentum 



components spatially which can then be measured by ab- 
sorption imaging 4\. This will project the sample into a 
definite component of the angular momentum and, due 
to the anti-correlation, the other sample will be projected 
into the opposite component. 

The low energy collisions of spin-1 atoms trapped in 
an optical potential can be described by the interaction 
potential: 
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where clf is the scattering length in the total spin F = 
0,2 channel and m is the atomic mass [12^. Here Pp is 
a projection operator which projects the pair of atoms, 
1 and 2, into a total hyperfine spin F state. Follow- 
ing the initial proposal of Fcdichev et al. 0], recent 
work has shown that the magnitude and the sign of the 
scattering leng ths can be tunned via all optical Feshbach 
resonances |l4l ITU Il6l Il7| . The enhanced experimental 
control this provides opens up a new arena of possibili- 
ties for quantum state control in these spin systems, of 
which, the present work is one example. 

We consider the case of a symmetric double-well po- 
tential (such as that recently realized at MIT jT||) where 
the system can be described in terms of well-localized 
spatial modes in the two wells of the potential each con- 
taining N/2 atoms. In addition, the spin dependence of 
the collisions is assumed to be small, \oq — CJ2I <C |ao|> 
so that we can assume the spatial modes are the same 
for each spin state [T2I Il9| . Defining the quantities 
Co = 4:Trh 2 (ao + 2<22)/3m and C2 = 4irfi 2 {a2 — ao)/Sm 
the Hamiltonian of the system has the form 
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where a,- iCe is the annihilation operator for an atom 
with spin projection a in the ith well. Here J — 

- J d 3 riPl(r)[-h 2 V 2 /2m + V pot {r)]Mr) and U F = 
Cf J d 3 r\ipi(r)\ 4 are given in terms of the spatial modes 
For each well we have defined the number op- 
erators hi = ^2 a a\ a ai t a and the spin vectors Lj — 

Lj <z z + \{Lj. + + L 3 -)x + - Lj^)y, where z, x, 

and y are unit vectors in a Cartesian coordinate system 
and 

Lj. z = aj^hj^i - a^^&j-x, (8) 
4- = V^tat _ x a ii0 + at Q a jA ) (9) 

and Lj,+ = Lj _. The above "tight-binding" Hamilto- 
nian is valid when the hopping can be treated as a weak 
perturbation to two independent wells, i.e. the hopping 
energy must be much smaller than the mode spacing in 
the independent wells |2fJ. 

When J = 0, it is convenient to introduce the simul- 
taneous eigenstates of hj, L 2 and Lj, z '- 

\ nj ,l 3 , mj ) ex (L i ,_)^-^(i])^(a] !l )^|vac), (10) 

where the spin-singlet creation operator is defined by 
A\ = a\ 2 a — 2aJ 1 aJ_ 1 |2l|. The singlet number Sj sat- 
isfies the relation 2s j — rij — lj . More generally, we can 
derive the operator identity 

hi(hi + l) =L 2 + A\Ai, (11) 

which has the interpretation: atoms either contribute to 
the magnitude of the spin angular momentum or they 
form singlets. It is also useful to define the simultaneous 
eigenstates of hj, L 2 , the total angular momentum L = 

L\+L\+2Li-L 2 andi z = Li^+L 2 , z - ^i, n 2 , Z 2 ; /, m), 
which are related to \n\, lx, m\)\\n 2 , I2, ^2)2 via the usual 
Clebsch-Gordan coefficients. 

In the limit of vanishing hopping (J = 0), we can 
use the fact that the L 2 commute with the Hamilto- 
nian to determine the ground state. In the case of anti- 
ferromagnetic interactions (C/ 2 > 0), the ground state 
of each well is the state that minimizes (L 2 ) [T2I I22] ] 
which, from Eq. (|llfl . implies that all atoms form singlets: 
(A\Ai) w (hi(hi + 1)). In the simplest case when N and 
ii\ = n 2 = N/2 are all even, the ground state is given by 
the product state 

|AF) = \h = 0, mi = 0)i|I 2 = 0, m 2 = 0) 2 (12) 

or |AF) = |Zi = 0,l 2 = 0;l = 0,m = 0), where 
we have suppressed the nt labels for notational con- 
venience. Using similar arguments we see that in the 
case of ferromagnetic interactions (U2 < 0), the angu- 
lar momentum in each well takes on its maximum value: 
(Lf) = (hi(hi + 1)). Note that in both the ferromag- 
netic and anti-ferromagnetic cases there is no correlation 
between the two wells. 



We are now in a position to describe the formation of 
the maximally entangled state |B/) using simple conserva- 
tion of angular momentum arguments. Assume our sys- 
tem has been prepared in the above anti-ferromagnetic 
(U 2 > 0) ground state such that I = 0. Now let us 
use an optically-induced Feshbach resonance to tune the 
scattering lengths ao and a 2 independently, such that U% 
is tunned adiabatically through zero until U 2 < 0. It 
is possible to do this far from resonance so the atom 
loss due to spontaneous emission from excited states can 
be minimized 1 15j (the experimental requirements of this 
transition are discussed in more detail below). Since the 
angular momentum is conserved in this process we have 
(L\ ■ L 2 ) = —{Lf) = —{Lf), but, now as we reduce the 
hopping coefficient, just as in the ferromagnetic case, the 
system must maximize the magnitude of spin in each 
well i.e. (Lf) = (n,-(nj + 1)). The only way to sat- 
isfy both requirements is for the spins to be maximally 
anti-correlated: (L± ■ L 2 ) = — {hi(hi + 1)). In the limit 
of vanishing hopping (for the simplest case of N even), 
m = ri2 — N/2 and the ground state is given by 

IBjv/s) = |Ii = N/2, l 2 = N/2, 1 = 0,m = 0), (13) 

which is the spin-Z singlet state. In terms of the an- 
gular momentum states of each well this state has the 
remarkable form given by Eq.J2J where we have sup- 
pressed the rii and U labels. The two wells are maxi- 
mally entangled as the reduced density matrix formed 
by performing a partial trace over one of the wells, 
p = ptr{|Bjv/ 2 )(Bjv/2|}, yields a completely mixed state 
with N + 1 degenerate eigenvalues 0. Generation of en- 
tangled states in two-level systems via adiabatic cross- 
ings similar to the present method have been discussed 
in Ref.||. 

In order to create an entangled state in a real dynami- 
cal process we must provide a coupling between the initial 
state |AF) and the final entangled state. Hopping be- 
tween the wells is one possible candidate to provide this 
coupling. Let us consider the hopping to be a weak per- 
turbation to two independent wells (i.e. J |J7 2 | <C Uo). 
To first order the hopping does not change the number in 
each well and the effective hopping Hamiltonian is given 
by 

Hgp = PH hop {E^ - flrito) _1 flhapA (14) 

where P is a projection operator onto the subspace of 
fixed number at each site and, since \U 2 \ < U , £7(°) m 
UqN(N/2 — l)/4 is the energy of the ground state with- 
out hopping. A similar approach to treat hopping was 
taken in Ref.0 for the case of two species in a double- well 
potential. In contrast to their scheme, where an entan- 
gled state is created from a highly excited state (different 
species in each wells) via hopping between the wells, our 
approach involves keeping the system close to the ground 
state and thus avoids decoherence due to excitations. Af- 
ter some straight-forward algebra (and dropping the con- 



3 



stant terms) we arrive at the effective Hamiltonian 
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this is the creation operator of a spin-1 Bell state) and 
U 2 = U 2 + 2J 2 /U . In the form given by Eq. (|16fl we 
have dropped the term proportional to the total angu- 
lar momentum as it is a constant and plays no role in 
the dynamics. Note that Ea. (|15l) shows that the hop- 
ping gives rise to a spin-spin coupling between the wells, 
which should be contrasted with a spin-spin coupling 
originating from a (usually weak) dipole-dipole coupling 
[24| . The energy states of interest are the eigenstates 
of -ffeff in the L 2 — 0, fij — N/2 subspace. When 
J = 0, for the symmetric case, the "bare" eigenvalues 
and corresponding eigenstates are E^ = — U' 2 lx{l\ + 1) 
and \h, l 2 = h,l = 0,m = 0), where h = 0,2,- •• ,N/2 
for N/2 even. These eigenstates are the generalized Bell 
states given in Eq. (J2J| as we can write 

|B (l ) - |Ji,Z a = ii,J = 0,m = 0) (17) 

h 

= {2h + l)- 1 ' 2 J2 (-i) m \h,m)i\h,-mW) 
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using the Clebsch-Gordan coefficients. For l\ > they 
are therefore maximally entangled states. A measure of 
the entanglement of a pure state is given by the entangle- 
ment entropy: S = — tr{plog 2 For the states | J B/ 1 ) 
we find that the entanglement entropy scales as the log 
ofh: S= log 2 (2Zi + 1). 

Using the effective Hamiltonian (|16f) we can analyze in 
detail the dynamical creation of an entangled state from 
the initial state |AF) in the case when U 2 {t) is varied in 
time as U 2 (t) = —at, where a > is a constant. For 
J = the energies of the Hamiltonian Ijltjfl undergo a 
multi- level crossing at U 2 (t = 0) = 0. For J 7^ the term 
proportional to &G couples these "bare" eigenstates and 
turns the level crossing into the avoided crossing shown 
in FigUl 

Numerical diagonalization (up to N = 120) of the 
effective Hamiltonian in the |B/ X ) basis shows that the 
minimum energy gap between the ground and first ex- 
cited state occurs at U 2 {t = 0) = and has the value 
A = 16J 2 /[/ for N > 1. Close to the adiabatic 
limit we expect that only the two lowest energy states 
play a role in the dynamics of the system and we can 
approximate the evolution as a Landau-Zener crossing 
|25|. Numerical simulations (see FigEJ) confirm that 
close to the adiabatic limit the probability of transitions 
to the first excited state follow the Landau-Zener law 
P = e' 2 ^, where T = A 2 / (hdE / dt) . Here dE/dt is 
proportional to the slope of the difference between the 
ground and first excited state energy levels and has the 
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FIG. 1: Example plot of the energy levels of the system as a 
function of t/^. The thin dashed lines depict the J = case 
and the solid lines the J 7^ case. The energy offset between 
the two sets of energy levels has been altered for clarity. 



form dE/dt = af(N), where f(N) is a function of N and 
satisfies f(N) < N(N/2 + 1) [c.f. FigHJ. The adiabatic 
limit is then given by a < 2A 2 /[KN(N + 2)]. Figure 
shows the evolution of (L 2 ){t) = —(Li ■ L 2 )(t) for dif- 
ferent values of a. From this figure we can see that for 
non-adiabatic changes the system tends to oscillate be- 
tween the different eigenstates of the Hamiltonian l|16fl . 

As m — mi + m 2 = 0, the states \Bi) are unaffected 
by a uniform magnetic field, however, a magnetic field 
difference between the wells leads to a term of the form 
fffieid ex Lx <z — E 2 z in the Hamiltonian which will break 
the rotational symmetry and allow (L 2 ) to vary in time. 
Spin flips and atom loss also give rise to changes in (L 2 ). 
In general, angular momentum conservation determines 
that a highly entangled state will form if the total spin 
angular momentum can be kept microscopic {(L 2 ) ~ 1) 
during the time that (L 2 ) is increased from ~ 1 to ~ N 2 . 

Once the entangled state is formed and the two wells 
are separated (J — > 0), applied magnetic fields will only 
lead to a phase rotation of the |Z,m)i|Z, — m) 2 states but 
will not effect the level of entanglement between the two 
wells. Two other important decoherence mechanisms for 
the entangled state |B;) are spin flips and atom loss. The 
state after a single spin flip in one of the wells is IV'flip) °^ 
L\ - IB;) which can be written as 



IV4 P }« E i-ircL-m+ilhrn-lUlh-mh (19) 
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where C l x y = + x)(l + y). Since the angular mo- 
mentum state of one well can be written as 
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FIG. 2: Growth of the spin angular momentum in well 
1 as U%{£) is varied from positive to negative from an ini- 
tial anti-ferromagnetic state. In this plot N = 120 and 
J = l?o/20. The solid line corresponds to 2ttA 2 /(M = 3200, 
the dashed line to 2n / (ha) = 320 and the dash-dotted line 
to 27rA 2 /(fta) = 32. The dotted lines show the angular mo- 
mentum for the three lowest energy eigenstates. The inset 
shows the overlap of the time evolved state with the excited 
state for the 27rA 2 /(te) = 3200 case. 



the annihilation operator acting on this state gives 



ai t i\l,l,m)i = 
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The state after the loss of one atom from the entangled 
state is \tpl oss ) K Si,i|Bj) which has the form 
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Since (B ; |V4 P ) = (B;|V>L S ) = 0, |B Z ) is not robust to 
spin flips or atom loss, as a single spin flip or the loss of 
a single atom will transform it into an orthogonal state. 
However Eq. (|19fl and Eq. (|24|) show that the state after a 
spin flip or the loss of one atom is still highly entangled, 
and therefore can still be a useful source of entanglement 
for quantum teleportationQ, for example. In figure|3]we 
show the reduction of entanglement entropy after a spin 
flip or atom loss as a function of the total initial number 
of atoms. These results suggest that weak decoherence 
has only a small effect on the entangled state. This is 
due to the fact that the decoherence effects on the level 
of the microscopic constituents (atoms) has little effect 
on the collective entanglement contained in the general- 
ized Bell state. These results also demonstrate that even 
without exact symmetry between the wells (slightly dif- 
ferent numbers of atoms in each well, for example), highly 
entangled states can still be created with this system. 

We now turn to a consideration of the general ex- 
perimental requirements for using an optically-induced 
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FIG. 3: Effects of decoherence due to the loss of one atom 
and one spin flip on the entanglement as a function of total 
number. The solid line shows the entanglement entropy of the 
state |Bjv/a), the dashed line that corresponding to IV'mp 2 ) ana - 
the dash-dotted line that corresponding to IV'ioss )• 



Feshbach resoance to tune the scattering lengths such 
that U2 a «2 - ao changes sign and the interactions 
change from ferromagnetic to anti-ferromagnetic or vice 
versa|l2|. Optically-induced Feshbach resonances work 
by off-resonant photoassociation of two free atoms to an 
excited molecular state, thereby altering the scattering 
length of the atoms [l3j. Since, in general, the excited 
molecular states with total spin F — and total spin 
F = 2 arc split by the hyperfine energy, hujhf = E2 — Eq , 
(see, for example, |2(|) it is in principle possible to tune 
the scattering lengths ao and ai independently as the 
photoassociative transitions to the excited molecular lev- 
els will experience different detunings. Following Ref.[l4| 
we can determine the change in the difference between 
the scattering lengths for a given detuning of the laser 
from resonance to be 
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where Ao is the detuning of the laser from the F = 
molecular state, E = h 2 k 2 /2m is the atomic kinetic en- 
ergy, 7 is the spontaneous emission rate and T is the rate 
of photoassociation and is proportional to the product of 
the laser intensity and the overlap integral between the 
free atoms and the excited molecular state. This equa- 
tion shows the "leverage" due to hyperfine splitting which 
enables us to tune the difference between the scattering 
lengths. 

Near an optical Feshbach resonance, spontaneous emis- 
sion of molecules in the excited molecular state inevitably 
leads to loss of atoms from the system. In the present 
case this loss rate is given by the sum of losses due to the 
two molecular levels and has the form 
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(24) 



For simplicity, in Ea. (|23|) and (|24|l . we have taken Aq ~3> 
Whf and set the rates 7 and T to be the same for the 
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two excited molecular states. In contrast to the experi- 
ment described in > nere we are n °t interested in large 
changes in the absolute value of the scattering length, but 
instead, in small relative changes in and ao, which can 
be achieved far from the resonance with a high laser in- 
tensity thereby keeping the loss rate to a manageable 
level To illustrate this let us consider the specific 
case of 7 Li (which is ferromagnetic with the scattering 
lengths ao = 12 a.u. and 0,2 = 5 a.u.[2^1 and a photoas- 
sociative transition from the atomic continuum of the 
ground molecular potential 3 S+ to the v — 72 vibra- 
tional level of the excited molecular potential 1 3 S+, as 
discussed in Ref. Q . The hyperfine splitting between the 
F = and F = 2 state in the excited molecular potential 
is LUhf — 274.5 MHz in this case [2^. For a laser intensity 
of / = 5 kW/cm 2 and temperature of lfj, K, Aa becomes 
positive (and the interactions become anti-ferromagnetic 
as required) for detunings A < 2.5 GHz to the red of 
the F = level. At this detuning, losses occur at the 
rate of i^i no i = 4x 10 -12 cm 3 /s, which is of the order of 
the rate of other loss mechanisms and is therefore of an 
acceptable level. 

In summary, we have shown that starting with two 
weakly-coupled anti-ferromagnetic Bose-Einstein con- 
densates it is possible to form a highly entangled state 



consisting of two anti-correlated ferromagnetic conden- 
sates by changing the sign of U2 via optically-induced 
Feshbach resonance, while conserving the total spin- 
angular-momentum. The entanglement becomes maxi- 
mal for symmetric wells in the limit of adiabatic changes. 
In addition to the theoretical significance of this entan- 
gled state (as a generalization of the Bell state) and its 
novel method of creation (via an anti-ferromagnetic to 
ferromagnetic transition) we have shown that the quan- 
tity of entanglement of this state is approximately pre- 
served under spin flips and loss of small numbers of atoms 
and therefore could be a useful resource for quantum in- 
formation applications. 

An interesting generalization of this idea is to an opti- 
cal lattice with ~ 1 atoms in each of the sites and hopping 
between neighboring sites|28|. Due to the number depen- 
dence of the adiabatic criteria described for the double- 
well case above, an entangled state in the lattice case 
may in fact be more easily realized experimentally. A 
detailed analysis of the dynamics of this system will be 
the subject of future work. 

The authors would like to thank F. Morikoshi and K. 
Shimizu for useful discussions and H. Pu for reading this 
manuscript. 
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